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Abstract 

Six-jet events via WW pairs, e+e“ —>• W~^W~ <liIhQ. 2 l !299 ci-re studied 

at tree level using helicity amplitudes. This is the dominant production 
mechanism for six-jet hnal states at Lep II energy. ISR effects are taken 
into account. Total production rates as a function of ycut are given. The 
relevance of these processes for the issue of colour reconnection is discussed. 
The cross section for hve-jet production via WW pairs at Lep II is also 
presented. 
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Introduction 


Several different mechanisms contribute to six-jet production. At Lep II energies they 
fall into two broad classes: the ‘point-like’ annihilation of to which domi¬ 

nates at LEP I energy, and vector boson pair production whose importance grows with 
increasing energy and dominates above the WW treshold. At still larger energy one 
should also consider the contribution of three vector bosons. Interference effects be¬ 
tween the two set of diagrams are expected to be small since the point-like contribution 
is dominated by the radiation of four gluons from the two primary quarks with only 
a small fraction of events with four and possibly six quarks. On the contrary in WW 
decays there are at least four quarks. Moreover only half of the WW events are flavour 
neutral and can therefore interfere with the point-like contribution. 

In this paper we study the production of off-shell WW pairs and their subsequent 
decay to hve and six jets at colliders. Since as{Mz) ~ .12, gluons are radiated 
with high probability and hve- and six-jet fractions are large also at experimentally 
relevant values of i/cut- This mechanism is the main source of hve- and six-jet events 
with total hadronic mass much greater than the mass of the Z. A large number of 
events with total hadronic mass close to Mz is produced by radiative return to the Z 
peak via initial state radiation. 

One of the main goals of Lep II is the measurement of the W mass with high accu¬ 
racy. Fully hadronic decays of the WW pairs represent 4/9 of the cross section. Apart 
from all the usual uncertainties related to measuring jet energies and directions, which 
can be partially eliminated by using energy-momentum conservation and possibly the 
approximate equality of the W~^ and W~ masses, there are uncertainties which stem 
from the fact that two decays occur in the same event. A jet from, say , the decay 
of the W~^ can be closer, in a given reconstruction scheme, to a jet from the decay of 
the W~ than to any other jet from the positively charged W. If the stray jet carries 
large energy, then the reconstructed masses will be quite far from the W mass which 
is already at present known with an error of about 200 MeV. One could imagine that 
such events could be discarded. However, due to the intrinsic width of the W and 
to all experimental uncertainties, it is impossible to impose very stringent cuts on the 
difference between the measured mass and the true mass. It is therefore important to 
estimate the effect of misassigned jets, including the large angle emission tails. For 
this purpose it is well known that matrix element calculations describe the full angular 
distribution of jets better than parton shower models while the latter are superior for 
small angle radiation. It should be mentioned that, due to the number and complexity 
of the diagrams involved, it is unlikely that a full 0{al) calculation will be produced 
in the foreseeable future, and only tree-level matrix elements will be available. 

Recently it has been pointed out i I, i that cross-talk between the decay of 
the two IF’s can take place also at energy scales much smaller than those typical of 
jets. This destroys the notion of two separate decays and cast doubts on our ability 
to reconstruct the masses of the original sources from the decay products, producing 
potentially large uncertainties. The perturbative contribution to the phenomenon, is 
related, in the Feynman diagram language, to the presence in the expression for the 
cross section of interference terms, in which real or virtual partons emitted in the decay 
of one W are reabsorbed by the decay products of the oppositely charged W. One of 


the simplest example is given by the interference between diagrams of type Ai and 
diagrams of type A 2 or between diagrams of type A 3 and diagrams of type A^ in £g. 1. 

The contribution of hard gluons, which can be treated perturbatively, to colour rear¬ 
rangement effects is expected to be small. The argument can be summarized as follows. 
Only interference terms between different diagrams can produce colour rearrangement; 
terms in which the gluons emitted by one line are reabsorbed by the same line do not 
mix the original decays. It is easily realized that all diagrams with a single gluon ex¬ 
change between the two fermion lines originating from different IT’s vanish when the 
sum over colours is performed. When two, real or virtual, gluons are exchanged, on the 
contrary non-zero contributions are possible. As a consequence colour rearranged final 
states are suppressed by compared with normal ones. Moreover, it can be shown 
that the color factor for the contributions where color rearrangement takes place is 
smaller by a factor {N^ — 1)“^ with respect to the color factor of the leading terms in 
1 /Nc at the same order in If one compares these terms with the lowest order cross 
section the ratio of colour factors is = |. A further suppression of the hard 

gluon contribution to colour reconnection derives from the hnite width of the IT’s. As 
an example consider the interference between the diagrams of set Ai and the diagrams 
of set A 2 in £g. 1. If the gluons carry a substantial amount of energy at least two of the 
W propagators must be far off mass shell and therefore the contribution of this term 
is negligible. The authors of ref. p]] give the following estimate for the ratio between 
the contribution of colour rearranged diagrams and the total cross section: 

a - \ N, J Mw ^ ’ 

In this expression is to be evaluated at = T^ where ~ .25 so that /Aaja < 
10-T 

The non-perturbative contribution to colour reconnection can be described in very 
simple terms in the string model, arguably the best tool for understanding the frag¬ 
mentation and hadronization phase of jet development, which links the perturbative 
gluon cascade to the observed hadrons. In this model, particle production results from 
the iterative splitting of color strings or antennae which are color singlets and evolve 
independently of each other. When two color singlet sources, like the two IT’s at Lep II, 
decay close to each other in space-time there two are different ways of forming colour 
singlets out of the two qq pairs. If the primary decays are IT’'' —and IT“ — q^q^ 
the qiq 2 and q^q^ pairs are always in a relative singlet state. However, with probability 
1/9, the qiql and q^q^ pairs can also form colour singlet states. Since the evolution of 
the antennae is determined by their invariant mass the multiplicity and distribution of 
detected particles depend strongly on the pairing pattern which is selected. 

A more quantitative estimate of color reconnection effects can be obtained using 
JETSET [^. In JETSET the primary decay is followed by a perturbative shower evo¬ 
lution, in which successive emissions are strongly ordered in angle as dictated by colour 
coherence. The showers develop as a Markov-chain process and no colour reconnection 
takes place in this phase. The two showers expand and when the hadronization phase 
begins they can overlap. During fragmentation each gluon acts as a qq pair and real 
hadrons are produced by the recursive splitting of colour strings extending between 
qq pairs. It is at this non-perturbative stage that the authors of ref. [Q allow colour 





strings to interact, to break and form new strings which may join quark originating 
from different showers, producing particles which cannot be assigned unambigously to 
any of the original hh’s. The details of the reconnection mechanism are unknown and 
require some additional modeling with respect to the well tested basic rules used in 
JETSET. It is this lack of information in the choice of the reconnection mechanism, 
that produce the spread in magnitude and even in sign of the predictions for the effects 
of colour reconnection in the measurement of My/. 


Results 


We have studied the reaction > W^W~ —*■ qiqiq 2 ^gg at tree level. Six quark 

hnal states have not been included since, at colliders, when a pair of gluons is 

replaced by a quark-antiquark pair, summed over all allowed flavour combinations, the 
cross section typically decreases by an order of magnitude. Quark masses have been 
neglected since the contribution of 6-quarks is severely suppressed by the smallness of 
the Vbc element of the CKM matrix. The relevant diagrams are shown schematically in 
£g. 1. For simplicity the lepton part of the diagrams is not drawn and a sum over the 
three possible structures which describe e+e“ ^ WW is understood. The set labelled 
A1 (A2) includes the twenty-four diagrams in which both gluons are emitted in the 
decay of the W~^ iW~)- Analogously, in A3 we include the twelve diagrams in which 
one gluon, let us call it g^ is emitted in the decay of the while the second, g^^ is 
emitted in the decay of the W~. In A4 the positions of the two gluons are interchanged. 

All matrix elements have been computed using the formalism presented in ref. 
with the help of a set of routines, called PHACT |p, which generate the building 
blocks of the helicity amplitudes semi-automatically. In our experience this method 
is faster than others which are commonly used 0.I.0- Further gains in speed can 
be obtained avoiding subroutine and function calls. In order to achieve this goal the 
routines in PHACT instead of computing the numerical values of the different terms, 
write the corresponding FORTRAN code, which can then be compiled and run. In the 
formalism of it is easy to save every sub-diagram and then to reuse it several times. 
With this procedure we have generated a rather large piece of code, which however 
runs quite fast, and therefeore can be used in high statistics Montecarlo runs. 

The amplitudes have been checked for QCD gauge invariance. We have used Mz = 
91.1 GeV, Tz = 2.5 GeV, Myy = 80.6 GeV, Pw = 2.06 GeV, sin^(6'w) = -23, irib = 5. 
GeV, ttem = 1/128 and = .115 in the numerical part of our work. 

Our main results are presented in fig. 2 through 4. Fig. 2a shows the total cross 
section as a function of i/cut in the JADE [|T3 scheme for the three energies which have 
been agreed on or are under consideration for Lep II, that is ^/s = 175, 190 and 205 
GeV. At such energies the typical value of ycut we have studied, namely ycut = 5 x 10“^ 
corresponds to a jet-jet invariant mass above 12 GeV. In Fig. 2a initial state radiation 
(ISR) has been neglected. We see that for a fixed value of ycut the cross section decreases 
with increasing energy, while, in this range, the WW cross section is increasing. This 
is due to the dehnition of ycut- 
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which means that at larger center of mass energy larger jet-jet invariant mass are 
required for an event to pass the cut. On the other hand the mass scale of the decays 
is obviously given by the W mass and does not change with s. Since small invariant 
masses are more likely for two jets from the same W that for jets from different W^s one 
indeed expects the fixed-i/cut cross section not to grow as quickly as the cross section 
for e+e- ^ WW. 

In Fig. 2b we study ISR effects for six-jet production at ^/s = 175 and = 205. 
The results at ^/s = 190 have been computed but are not shown for the sake of clarity. 
We have used the structure function approach at leading-log as in ref. [O]. The cross 


section is decreased by about 20% at ^/s = 175 and by about 10% at -y/s = 205. The 
effect at ^/s = 190 nicely interpolates those at the two energies which are shown. 

For comparison in Fig. 3 we present the cross sections for five jet production at 
the same energies used for six jets. We have repeated the calculation of ref. []^| and 
found complete agreement. The results without ISR are given in fig. 3a and those 
which compare results with and without ISR effects in fig. 2b. The dot-dashed line 
in fig. 3a shows the five jet background cross section [^, Q from e+e“ —> qqggg and 
e+e" ^ qiWq 2 ^g- 

There is an aspect which deserves some comment. When our results for five-jet 
production [Q at Lep I are compared with the data presented by ALEPH and 
OPAL llTBl it is clear that the absolute normalization is about a factor of five too small. 


The simplest explanation for this discrepancy is our choice for In fact we have 
used as = .115 which corresponds to = M|o with A;^ = 200 MeV with five active 
flavours in the standard formula: 


as{Q) = 
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6olog(Q7A2) 
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hi log(log(QVA^)) 

hi log(QVA2) 
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The analysis of shape variables and jet rates to Oi^al) has shown that, in order to get 
agreement between the data and the theoretical predictions, the scale of the strong 
coupling constant has to be chosen to be Q = x^Mzo, with 0.1 [|^. It has later 


been shown that when the relevant logarithms are properly resummed agreement is 
obtained for much larger values of the scale, ~ 1. 0. It is therefore not surprising 
that our tree level expressions require a relatively small scale in order to describe the 
data. As a consequence the cross sections presented in Fig. 2a,b are expected to 
be somewhat underestimated. Since ag is an overall factor our results can be easily 
modified if a different value for the strong coupling constant is preferred. 

Leaving aside the issue of best choice of as, with an expected luminosity between 
300 and 500 ph~^ it is clear that only for i/cut < 1-5 x 10 "^ six jet events can be observed. 
The cross section grows very rapidly and for i/cut = 5 x 10 "^ one expects 0 ( 100 ) decays 
of a WW pair to six jets per year and per experiment. 

In Fig. 4 we present the gluon energy spectrum. The continous line gives the differ¬ 
ential distribution for one fixed gluon, while the dashed line and the short-dash-long- 
dash line give the energy spectrum for the most and least energetic gluon respectively. 
These quantities are not directly observable but since in most cases the softest jets 
are the gluon ones they give an indication of what sort of low energy jets one might 
expect in WW decays. The softest gluon energy distribution is rather narrow and 




























peaks at about 8 GeV. The most energetic gluon has a much broader distribution with 
a maximum at about 20 GeV. 

In order to obtain a quantitative estimate of the colour reconnection contribution 
to the cross section we have integrated over phase space the interference between the 
diagrams in the set Al and those in A2 and the interference between the diagrams 
in the set A3 and those in A4. The interference between (A1+A2) and (A3+A4) is 
zero. In the following we will call the result aint- These integrals are particularly 
challenging. The magnitude and phase of the integrand change rapidly at all W poles 
and since different sets of particle reconstruct the W mass in different diagrams it 
is impossible to eliminate the strong peaking structure with the standard change of 
variable kf — M^ = MyyTwtanOi. In addition the requirement that all jet-jet invariant 
masses be larger than = si/cut introduces discontinuities within the integration 

region. The accuracy of our results for the interference terms is typically of the order 
of l-i-2%, the least accurate point being at i/cut = 1 x 10“^ and = 2.06 GeV where 
the error is about 6.6%. The results presented in £g. 2 through 4 have an accuracy 
well below 1%. 

In the narrow width limit the cross section scales as (T^y)”^ and the bound (0) 
suggests that the interference terms are proportional to (Ty/)”^. Therefore we have 
investigated the behaviour of aint as a function of the W width in the range between 
10 and 2 GeV. Numerical instabilities become less severe if one artihcially increases the 
IV-boson width and this allows us to keep the numerical accuracy of the integration 
under control. Our result are presented in Fig. 5 for a/s = 175 GeV and two ycut 
values, i/cut = 6.5 x 10“^ and i/cut = 1 x 10“^. For comparison we also plot the 
curves of the form y = a/Vw + h which interpolate our results at Fyy = 10 GeV and 
Fvy = 7 GeV. In events which do contribute to atnt the gluon energies are restricted 
from below by the cut on invariant masses and from above by hF-width effects as 
explained previously. The resulting band becomes narrower for increasing ycut and for 
decreasing Tw This can be clearly seen in Fig. 5. For large values of the IF-width 
the results approximately behave as a/Vw + h while for values of F^y approaching the 
physical value they gradually fall below the reference curve. This is more evident for 
Ucut = 6.5 X 10“^ where the lower limit of the allowed band is higher. Taking the cross 
section for —>• W~^W~ to be 8 pb at a/s = 175 GeV, /S.a/a turns 

out to be 1 X 10“® at ycut = 6.5 x 10“^ and 5 x 10“^ at ycut = 1 x 10“^, within the 
bound of ref. |^. The ratio of the interference terms with the six-jet cross section has 
a milder dependence on ycuu being 3.3 x 10“^ at ycut = 6.5 x 10“^ and 5.7 x 10“^ at 
ycut lx 10 

Conclusions 

We have computed at tree level the cross section for the process e'''e“ —>■ W~^W~ —>■ 
qi(hq 2 ^gg which is the dominant contribution to six-jet production at Lep II energies. 
With a luminosity of 500 pb~^ one expects (9(100) decays of a WW pair to six jets 
per year and per experiment for pcut = 5 x 10“^. We have studied the lowest order 
non-trivial perturbative contribution to colour reconnection. In the range we have 
considered, namely y > 1 x 10“^, it is small and within the bound of ref. |^. We have 
presented the cross sections for the process e+e“ —> W^W~ qiWq 2 ^g at Lep II. 
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Figure Captions 

Fig. 1 Representative diagrams of the decay part contributing to e’''e“ — > WW —> 
QiQ 2 QsQ 4 : 99 - For simplicity the lepton part of the diagrams is not drawn and a sum 
over the three possible structures which describe e+e“ — WW is understood. In 
set A1 (A2) we include the eight diagrams in which both gluons are emitted in 
the decay of the W^ iW~). Analogously in set A3 we include the four diagrams 
in which gluon is emitted in the decay of the W^ while is emitted in the 
decay of the W~ . In A4 the positions of the two gluons are interchanged. 

Fig. 2 Total cross section as a function of Hcut in the JADE scheme for six-jet pro¬ 
duction via WW at = 175 GeV (continuous line), ^/s = 190 GeV (dashed 
line) and -/s = 205 GeV (dotted line). In Fig. 2a ISR effects are not included. 
In Fig. 2b we compare results which include ISR effects with those without ISR, 
as explained in the main text, at = 175 GeV and = 205 GeV. 

Fig. 3 Total cross section as a function of ycut in the JADE scheme for hve-jet produc¬ 
tion via WW at ^/s = 175 GeV (continuous line), ^/s = 190 GeV (dashed line) 
and = 205 GeV (dotted line). The dot-dashed line in Fig. 3a gives the QGD 
contribution to hve-jet production through point-like annihilation at ^/s = 175 
GeV. In Fig. 3a ISR effects are not included. In Fig. 3b we we compare results 
which include ISR effects with those without ISR, as explained in the main text, 
at ^/s = 175 GeV and ^/s = 205 GeV for the WW signal only. 

Fig. 4 Gluon spectra in six-jet production via WW at y/s = 175 GeV. The conti- 
nous line gives the energy distribution for a single gluon, while the short-dashed 
line and the long-dash-short-dash line give the spectrun of the most and least 
energetic of the two gluons, respectively. 

Fig. 5 Integral over phase-space of the sum of the interference between the diagrams 
in set A1 with those in set A2 and of the diagrams in set A3 with those in set A4 
at ^/s = 175 GeV. 









6 jet distribution 



Fig. 2b 










5 jet distribution 



Fig. 3a 






(qd) (^)D 


5 jet distribution 



10 '* 


No ISR 








dff/dE (pb/GeV) 


6 jet: Gluon en 


ergy distribution 



Fig. 4 








